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A method is described for predicting the fast fading rate of moon-reflected signals. It is 
based entirely upon considerations of the observer-moon positions and relative motions. 
Experimental results which are in good agreement with the computed fading rates have been 
obtained from a moon-reflection experiment at a frequency of 412 megacycles per- second. 
Some possible implications of this method of interpreting fading rates are given. 



1. Introduction 

The effect on UHF signals of transmission through 
the ionosphere and into space is becoming a problem 
of increasing importance 1 . At present the 4 char- 
acteristics of the transmission medium are not known 
in detail, and the effects upon transmitted signals 
remain questionable. 3 It is known that under certain 
conditions, e.g., auroral activity, the disturbances 
are such that radar returns may be obtained from 
them. While this provides valuable information 
about the backscatter which may be expected, it docs 
not give a direct measure of the effect upon the trans- 
mitted signal. Some information may be obtained 
by making use of the signals received from radio 
stars; the scintillation rates of such signals have been 
observed to increase at times of disturbed ionospheric 
conditions. However, a more controllable extra- 
terrestrial signal source has some advantages com- 
pared with radio stars. Prior to the advent of 
artificial satellites, it appeared that use of the moon 
as a reflector was a reasonable means of obtaining 
such a signal. Moon-reflected signals have been 
reported for some years [2], 4 and the use of moon 
relay circuits for communications purposes is de- 
veloping, but most of the measurements reported 
have not been directly applicable to the problem of 
ionospheric transmission effects upon UHF signals. 
The results reported in this paper were obtained in 
the course of a moon-reflection experiment carried 
out between South Dartmouth, Mass., and Alpha, 
Md. [3]. The procedures established in this test 
formed the basis for a later transmission experiment 
carried out between College, Alaska, and Westford, 
Mass. It is hoped to publish a full account of this 
latter experiment at a later date. 

Among the results of the test was the direct 
correlation between the fading rate of the received 
signal and the predicted calculations of the total 
libration rate of the moon. The spectral broadening 
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of the lunar signal was thus shown to he an effect of 
the motion of the moon. The major portion of this 
paper is concerned with an outline of these calcula- 
tions and a demonstration of their agreement with 
measured values. 

2. System Used in the Experiment 

The transmitter and 60-ft diam parabolic antenna 
were located at the M.I.T. field station at 'Round 
Hill, South Dartmouth, Mass. (41.5395° N, 70.9512° 
W). The receiving site was at Alpha, Md. (39.3224° 
N, 76.9258° W), about 350 miles from South Dart- 
mouth. The receiving system used a 28-ft diam 
parabolic antenna with two orthogonal feeds. Two 
identical receiving systems were used with common 
local oscillators. With a transmitted power of 40 k v 
at 412 Mc/s and a final receiver bandwidth of 50 c/s, 
a signal-to-noise ratio of 25 to 30 db was obtained. 
Paper tape records were made of the detected signals, 
and periodic magnetic tape recordings were made of 
the signals at an intermediate frequency of 2,500 c/s. 
Observations were made from moonrise to moonset 
each day from 6 August until 29 August 1957. 

3. Paper-Tape Recording for Fading-Rate 
Measurements 

Typical recordings of the detected signals are 
shown in figures 1, 2, and 3. The frequency of the 
"fading" shown is a prominent feature, and during 




Fig 



ure 1. Fast recording of detected signal, slow fading rate. 
(Traced from original recording chart.) 
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Figure 2. Fast recording of detected signal, medium fading rate. 
(Traced from original recording chart.) 




Figure 3. Fast recording of detected signal, fast fading rate. 
(Traced from original recording chart.) 



the test it was found to vary within fairly wide 
limits. At times the rate appeared to be several 
cycles per second, at other times the signal level 
remained almost constant for several minutes. 
Figure 1 shows a typical section of a recording taken 
at a relatively slow fading rate period, figure 2 shows 
an average fading rate, and figure 3 shows one of the 
faster fading rates. All of these were taken during 
the test with a detector time constant of 0.15 sec. 
The magnetic-tape recordings of the 2,500 c/s IF 
signals formed a convenient stored source of signals 
to be played back through a detector with an 
adjustable time constant. Some experimentation 
was needed to determine a suitable time constant 
which, it was felt, produced a record that showed up 
most of the signal variations without being too 
affected by noise. A value of 0.015 sec for the time 
constant was chosen as a suitable compromise, and 
many of the recorded 2,500 c/s IF signals were 
replayed with this modified detector system. Typical 
results obtained on a Brush recorder are shown in 
figures 4 (a) and (b). The paper-tape recordings 
thus obtained were used to measure the fading rate in 
a rather arbitrary but consistent manner. This was 
done simply by counting the number of detected 
signal maximums occurring over an interval of 
approximately 1 min at any given period. "Results 
varied from 3 to 4 fades/sec to approximately 0.005 
fades/sec. The higher frequencies were observed 
most of the time, while the very low fading rates 



were comparatively infrequent, and lasted for only a 
short time in any case. Some of the measurements 
are shown plotted as fading rate versus time of day 
in figures 11 (a) to (g). 

4. Libration-Rate Computations 

It seemed to be a reasonable assumption that the 
effective libration of the moon had a large part to 
play in producing the varying fading rates of the 
moon-reflected signals [2]. (The computations to be 
described first determine what is termed here the 
"total libration rate" of the moon, and thus proceed 
to calculate the resulting "Doppler spread." This 
can then be interpreted as a fading rate, and compared 
with the measured fading rate.) 

An overall view of the geometry involved in the 
calculation may be useful. Figure 5 shows a simpli- 
fied sketch of the earth-moon system. In effect, the 
moon presents a slightly varying face to the earth, 
due to the dynamics of its motion. This means that 
an observer on earth looking toward the center of 
the moon does not always see the same surface point 
(A in fig. 5), but rather a changing point on the surface. 
The motion involved is described by the selenographic 
latitude and longitude figures listed in the American 
Ephemeris and Nautical Almanac. The axes about 
which these positions are measured are shown in 
figure 5. In general, the axis for the selenographic 
latitude variations is inclined to the celestial equator. 
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Figure 4 a, b. Fast recordings of detected signals with 0.015-sec, time constant. 
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Figure 5. Geometry for fading rate computations. 
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An observer on earth, at point B, must look in the 
direction BA in order to see the moon. Then, due 
to the earth's rotation about the axis OP, the ray BA 
sweeps across the moon's surface. On addition, the 
combination of these two effects give rise to what is 
termed the "total libration" of the moon. A system 
of axes can be set up, and the motion of the moon 
can be resolved along the same set of axes. Conibi- 
nation of the two sets of components is then a simple 
1 natter. In effect, expressions are derived for the 
total libration rates in latitude and longitude, l T p and 
l TL , respectively, of the form 



and 






where V$ and l' L are libration rates due to the moon's 
motion, and/i and/ 2 are equivalent functions due to 
the earth's motion. All the functions, of course, are 
time-dependent. The total libration rate is given by 
the square root of the sums of the squares of 1 T b and 
l TL . Hence, both of these quantities must become 
small in order to give a low value for the total effective 
libration rate. 

With the libration rate computed, the next step is 
to transform it into an effective fading rate. Suppose 
that instantaneously the moon had only one strongly 
reflecting region located centrally at point A, so that 
a cw signal reflected from this region would show only 
the main Doppler shift in frequency. If now two 
additional reflecting portions are assumed to exist at 
C and D, toward the rim of the moon and symmetri- 
cally located with respect to the center portion, then 
each of these regions would be in motion due to the 
effective libration of the moon. The component of 
this motion resolved along the direction to the ob- 
server would give rise to a small additional "Doppler 
spread," so that an incident cw signal would be 
reflected with two Doppler sidebands. Reception 
and detection of this signal would then show a signal 
of varying level, with the Doppler spread in effect 
producing the "fading rate." Now actually a large 
portion of the moon must be effective in giving rise 
to appreciable reflections, and each small contributing 
area gives rise to its own particular Doppler spread. 
Integration over the whole surface is difficult, since 
the reflection mechanism is by no means well under- 
stood, and the moon is a rough body. However, it 
is possible to examine the moon's surface to determine 
the maximum extent of the Doppler spread. Con- 
sider a circle drawn on the moon's surface, the radius 
of this circle being specified as a given fraction of the 
actual full radius of the moon All points on this 
circle can then be examined to determine that point 
which gives rise to the maximum magnitude of 
Doppler spread. If the point C on the moon's surface 
is considered to be the point in question, then its 
position on the circle is specified by the angle v , meas- 
ured from the y' axis. It is apparent that major 
emphasis in the calculation must be given to a careful 
consideration of the geometry involved. 



4.1. Moon's Position 

Figure 6 shows a sketch of the celestial sphere 
with the celestial equator in the same plane as the 
earth's equator. The ecliptic has its ascending node 
on the celestial equator at the point y, the vernal 
equinox, while the ascending node of the moon's 
orbit on the celestial equator is marked by the point 
y\ The moon is assumed to be at the point B, so 
that the right ascension of the moon, BA m , is given 
by the arc yC, and the declination, 8 m , by CB. 
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Figure 6. Moon-position geometry. 





The angle 6 is measured in the plane containing 
the radii Oy' and OB, and is taken to be positive in 
the direction from Qy' to OB. The angle 8 t is the 
inclination angle of the great circle, centered at 0, 
and passing through the points y' and B. 

At any particular time the position of the point y' 
is given by noting the value of the moon's right 
ascension during the preceding month when its 
declination is zero, and when it is changing from 
negative to positive values Hence, with this value 
of right ascension BA in degrees, y'y is given by 



and 



y'y=360°-iLlo, 

y'C=y'y+BA m 



(2) 



Thus, from the triangle y' CB, the angle 6 is speci- 
fied by 

cos0 = cos(y'y+/?.4 m )cos(> m , (3) 

and 



8 t =2 arc sin 
where 



/ sin (S-0) sin (S -y'y- RAJ 
V sin 6 sin (y'y +RA m ) 



(4) 



S=[6 + y'y+BA m +\8 m \]/2. 

4.2. Orientation of Moon in its Orbit 

Figure 7 shows the moon's equatorial plane, with 
its ascending node on the celestial equator desig- 
nated as y"'. The angle 0' is defined as the distance 
along the celestial equator from the true equinox to 
the ascending node of the moon's mean equator, 
so that the arc y'y"' is given by 



7 7 



^' + 7'7 



(5) 



Consider two vectors A h A 2 , where the vector A { 
is parallel to the equatorial plane of the moon and 
perpendicular to OB, while A 2 is orthogonal to A t 
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Figure 7. Moon orientation geometry. 

and normal to the equatorial plane and directed 
upward. These vectors describe the axes about 
which the moon's selenographic latitude and longi- 
tude are measured. For computational purposes, 
it is convenient to transform these two vectors into 
two vectors A[, A 2 , where A[ is in the moon's 
orbital plane, perpendicular to OB, and A 2 is 
orthogonal to A\ and OB, Let three unit vectors 
be defined as follows: 

V=unit vector in direction OB, 

V=unit vector in direction OZ (OZ perpen- 
dicular to orbital plane), 

i p '=xmit vector in direction OP (OP perpen- 
dicular to equatorial plane). 

Consider an xyz-8Lxis system, with the z-axis coinci- 
dent with 01 and with the //-axis in the orbital plane. 
Then the angle e is defined as the angle from the 
x-axis to OB, measured in the orbital plane. The 
angle U is the included angle y f ly ,n (see fig. 7). 

Unit vectors i a , } a yf , in the direction of A\ and A[ 
respectively, are given by 



and 



"V /\ ' x' 



ly' l Z ' /\ fj'. 



(6) 



(7) 



The vectors A u A 2 can be described with the aid of 
unit vectors as 



and 



A 1 =A l i a 



A 2 =A 2 (—i y sin Zf-f-v c° s h)- 



(8) 



(9) 



The relations between the two sets of vectors A u A 2 
and A[, A 2 give the following expressions for the 
magnitudes A[ and A 2 : 



A[=(A l +A 2 )'i 9 



A, 



Vl + cos 2 e tan 2 l t 



-A 2 sin l t cos e 
(10) 



if / . i ~t \ -^i] cos € van L 7 . A 7 /. M v 

A 2 =(A 1 +A 2 ) ■ v= i~ 9 * +A 2 cos l t (11) 

yl + cos^ e tarn /* 

The angle i is defined as the inclination of the 
moon's mean equator to the earth's true equator. 
With knowledge of y'y' n t i and 8 t , the triangle 
y'y'"l can be solved for l t to give 

Zi=arc cos [cos h { cos i+sin b L sin i cos (y'y ,n )\ 

(12) 



and 



y'Z=arc cos 



tcos hi cos li— cos i~\ 
sin <5 Z - sin Z,- J 



(13) 



With the angle 6 obtained from (3) and y 7 fro in (13), 
the angle e is given by 



e=6-y'l. 
4.3. Effect of Nonspherical Earth 



(14) 



If a meridiana! cross section of the cart h is assumed 
to be elliptical in shape, then each point on the sur- 
face has a geocentric latitude that is slightly different 
from its normal geographic latitude. In addition, 
of course, the radius from the center point of the 
earth to the observer varies with latitude. The 
relations among the following parameters: 

a eQ = equatorial radius of earth, 

r= radius to point at geocentric latitude <£ c , 

6=ellipticity, 
c = geocentric latitude, 

0=geographic latitude, 



(15) 
(16) 



are given by 

r2== ael(l-e 2 ) 
1 — e 2 cos 2 <p e ' 

tan (p e =(l — e 2 ) tan <p, 

e 2 =0.006768658. 

4.4. Doppler Shift of Moon-Reflected Signal 



The sketch shown in figure 8 demonstrates the 
geometry involved. In essence, for an observer at 
position A, it is necessary to calculate the rate of 
change of the distance D. The moon's position is 
given by the distance 1) () and tin 4 coordinates shown 
in figure 9. It is necessary to note the distinction 
between the two sets of coordinates shown in this 
figure. The unprimed set is related to the earth's 
equatorial plane, while the primed system is related 
to the moon's orbital plane. The following param- 
eters are shown in figures 8 and 9. 
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D = distance from center of earth to center of 
moon, 

Z?=distance from observer at A to moon's 
center, 

a a = distance between earth's center and ob- 
server at A, 

\f/= angle between D Q and a a , 

5*=angle between earth's equatorial plane and 
the moon's orbital plane, 

y'C= right ascension of moon +7'y ; 

d m = declination of moon, 

<f) a = geocentric latitude of observer, 
y'D= right ascension of observer +t'y. 




Figure 8. Geometry for distance from observer to moon. 




MOON'S POSITION 

Figure 9. Geometry for positions of observer and moon. 

In order to express the angle \p in terms of the 
position parameters two unit vectors are introduced. 
A unit vector i a is defined in the direction OA, while 
ix is defined as a unit vector in the direction of D in 
figure 8 (OB in fig. 9). These vectors may be 
expressed as 

i a =i z sin (t> a +ix cos 4> a cos (y'D) + i y cos <j> a sin(y'D), 

(17) 
and 

ix =i z sin 8 m + i x cos8 m cos (y'C)+i y cos 8 m sin (y'C) 

(18) 



The angle \f/ is given by 
cos\l/=i x ' -i a 
.'. cos^=sin 8 m sin a +cos<£ a cos 8 m cos(LHA) 

where LHA is the local hour angle of the moon, and 
may be expressed as 

LHA= (y'D) - (y'C) = (yD) - ( y Q 
Now the rate of change of the distance D is given by, 

D=~^Dl+a 2 a -2a a D cosrP 



D- 



D (D —a a cos \f)—D^a a -^ (cos \f) ( . 



D 






The values of horizontal parallax, n, listed in the 
American Ephemeris may be used to compute D , 
giving 



Po= ihTn' 



and 



D =- 



a ea cos IT dU 



sin 2 n dt 



It is also necessary to differentiate cos \j/ from eq (19), 
which in turn involves the rate of change of the 
local hour angle. In this manner the angular rota-^ 
tion rate of the earth and the rate of change of right 
ascension of the moon are brought into the computa- 
tions. Thus 

-77 (cos \p) = 5 OT [cos 8 m sin a — sin 8 m cos <£ a cos (LHA)] 



and 



-cos 4> a cos 8 m sin (LHA) -y- (LHA) 



^LHA)=. e - Tt (RA m ), 



where co e is the angular velocity of the earth, and 
d/dt (RA m ) is the major component of the angular 
velocity of the moon about the earth. 

The above expressions enable D a and D b to be 
derived for the transmitting and receiving stations, 
so that the total Doppler shift i d is given by 



U=~k (4+ A) 



(21)1 
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where D a and D b now represent the distances between 
the transmitter and the moon and the receiver and 
the moon, respectively. 



4.5. Libration-Rate Computation 

The essence of this computation lies in the descrip- 
tion of the motion of a point on the moon's surface 
as seen by an observer on earth. Part of this motion 
is due to i he actual libration of the moon, and part is 
due to the movement of the observer. Thus the 
computations involve the rates of change of the 
selenographic latitude and longitude, plus the rate 
of change of some quantity which describes the 
geometrical earth-moon relationship. The effects of 
these two plrysically distinct processes, when added 
appropriately, then give the total libration rate. It 
is convenient to use the sign convention used in the 
American Ephemeris. Thus, the selenographic longi- 
tudes are measured in the plane of the moon's equa- 
tor, positive toward the west, the axis of reference 
being the radius of the moon which passes through 
the mean center of the visible disk. Latitudes are 
measured from the moon's equator, positive toward 
the north. 

In section 4.2, two vectors A u A 2 were considered. 
Now, let A\ become the axial vector ISpjind A 2 be- 
come the axial vector Sl, where Sp and S L represent 
axes about which the selenographic latitudes and 
longitudes are measured. The time derivatives of 
the scalar values of these quantities are given by 



ls= 



dt 



and 






h=- 



Also, in sections 4.2 and 4.4 primed quantities were 
used with reference to the moon's orbital plane. The 
same significance again is attached to primed quanti- 
ties here. Equations (10) and (11) can be rewritten 
to give Sp and S L , and the resulting expressions 
differentiated to give 



h 



dt ~ p Vl + cos 2 e tan 2 /, 



■ l L sin li cos e 



L(i 

la 



cos 2 € tan 2 /, sec 2 / ?: , ~ T ~]dli 

+ cos*«tan*/,)«* +& C ° S € C ° S U \ Tt 



'p cos e tm\ 2 l t sin e 



'l-fcos 2 e tam/,) 3/2 



' ml *]w 



(22) 



and 
dSL_ 7 , 

~TT — L L 



1$ cos e tan /,: 



Vl+cos 2 € tan 2 /, 
Sfi cos e sec 2 /* 



-f l L cos l t 



+[■ 



-[a 



S L sml^j t 

I* 

+ cos 2 e tan 2 / Z ) 3/2 J dt 



(1 + cos 2 e tan 2 /,) 3 / 2 

Ss sin e tan / 



(23) 



The effect of the observer's motion may be intro- 
duced through use of the parallax angle X}, shown in 
figure 8. Consider ^ to be an axial vector repre- 
senting a small rotation when the magnitude of X) 
represents the angle between D and D . Hence for 
small S> S is given by 



where 



2_j — \l x r • %a) jj 



Di=D —a a cos ^ 



(24) 
(25) 



In general, the observer's motion will give contribu- 
tions to both the longitudinal and the latitudinal 
libration rates. Bearing in mind the association of 
the primed coordinates with the moon's orbit, and 
the fact that the selenographic latitudes and long- 
itudes were adjusted to this set_of coordinates, it is 
apparent that the projections of ^ on the y' and z' 
axes are required. The component in the z' direction 
is given by 



(26) 



(27) 



(28) 



Z-jL — ( V X V) * % z' TT — jT ^a ' V 

and in the y f direction by 

From figure 9, it is seen that 

i 2 r = i z cos hi—i y sin 5, 
and the unit vector i y , is given b}' 

tyf ^=% z r /\ % x i 



with v and i x , described by eqs ^28) and (18) respec- 
tively. The dot products in eqs (26) and (27) may 
now be expanded to give the components of Xj as > 

XU^rf {cos hi cos d m cos 4> a sin (LHA) 

— sin 5, [sin<5 m cos</> a cosY'J9— cos 5 m cos (y'C) sin</>J} 

(29) 
and 

X^lr [sin 5 t cos a sin (y'D) — cos 5, sin a ] (30) 

The above two equations can be differentiated to 
give their rates of change. This involves a number 
of terms, as most of the parameters are time de- 
pendent. However y'y and 8 t change so slowly that 
negligible error is introduced by assuming them to 
be constant. The terms thus involve the rate of 
change of the distance, A; the angular rotation rate 
of the earth co e ; and the rates of change of the moon's 
right ascension and declination. Hence ^ L and y^ g 
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may be expressed as 

^-iL=—jy 2 L +u c jj cos <p a [cos 8 t cos 8 m cos (LH A) +&in 5 t sin 8 m sin (y'D)] 

—-7- (RA m ) cos <5 m y|- [cos 5 X cos <p a cos (LiL4)+sin 5,- sin <^ a sin (7'C)] 



~y7 5 m {sin <5 W cos <^ cos <p a sin (LHA) 
D 1 



-sin 5^ [cos <5 W cos <p a cos (7'/))+ sin <5 W cos (7 '(7) sin <p a ] 



(31) 



and 



A 



2*=-g Z*+«. ^ sin 5 2 cos ^ cos ( 7 'Z?), (32) 

where 

Di=D -a a -^ (cos$. 

In obtaining the sum of the effects of the moon's 
motion and the earth's motion, it must be remem- 
bered that a rotation in one sense on earth appears 
as the opposite sense on the moon. That is, an 
increase of the parallax angle toward the west or 
north appears as a rotation of the moon toward the 
east or south. Hence, the components of the total 
libration rate for a two-station system are 



/ , 1 • 1 • "> 

l -U'- l ± -is "V 

Itl — ^ II o La 9 I 



(33) 
(34) 



where the subscripts a and b refer to the transmitting 
and receiving stations. The total libration rate, l T , 
is thus given by 



/, 



~-'\It0-{-1'T1 



(35) 



4.6. Doppler Spread of Moon-Reflected Signal 

So far, when considering the geometry of the 
earth-moon system, the reference point on the moon 
has been its center. The corresponding point seen 
on the surface of the moon may be imagined to be 
the main reflecting point. However, as mentioned 
previously, many other points on the moon's surface 
are effective as reflectors. Any one such point will 
be in motion due to the libration of the moon, and in 
consequence of the fact that its line-of -sight velocity 
is different from that of the central point, its Doppler 
frequency will be slightly shifted with respect to that 
of the central point, i.e., with respect to the main 
computed Doppler frequency. This shift in Doppler 
frequency is referred to as the "Doppler spread." 
Obviously, this will depend upon the position of the 
offset region as well as on the angular velocities in- 
volved. Consequently, in order to obtain an idea of 
the maximum Doppler spread that may occur over 



a given portion of the moon's surface, it is convenient 
to specify a ring whose radius is a given fraction of the 
radius of the moon, and to examine this ring from 0° 
to 360° in order to pick out the maximum value of 
the Doppler spread. 

This method of computing the Doppler spread 
allows a number of approximations to be made. The 
choice of which approximations to make can be 
based upon the magnitude of the total libration rate. 
When this is not too small, the Doppler spread is of 
the order of cycles pers second, and any refinement 
in the computation would have a relatively small 
effect. As the libration rate decreases, so does the 
magnitude of the Doppler spread, and at some point 
the errors introduced by the approximate computa- 
tions may become significant. The overall accuracy 
depends, of course, upon the accuracy of the input 
information and, if this is not very precise, any small 
differences carried in the computations may not 
show up. 

Figure 10 shows a sketch of the geometry involved. 
The moon's center is at and the observer at 0'. 
Although only one observer's position is shown, the 
computations are carried out for two separated sta- 
tions, i.e., transmitter at station a and receiver at 
station b. The axes shown at the moon's position 
refer to the moon's orbit, as mentioned previously, 
and the quantities l TL and l T p are used to describe 
the moon's apparent rotation about these axes. The 
circle to be examined is shown as GLCH with radius 
FL, or kR m , where R m is the moon's radius and k is 
the fractional moon's radius. The position of the 
point L is specified by giving the value of k and the 
angle v, where v is measured positively upwards 
from a line parallel to the axis of l T p. 




Figure 10. Geometry for Doppler-spread computations. 
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Due to the apparent rotation of the moon about 
the l T p axis, the point L has a component of velocity 
directed towards 0'. The approximation is made 
that station a may be assumed to be al 0', so that an 
observer sees the incremental velocity of point /, as 

v a p=kR m smv(lp—2fia). 
Rotation about the l TL axis gives a similar component 

VaL=—kR m cosj> (V L — 2 La )> 

The Doppler spread for each station is due to 
these incremental velocities. The Doppler spread 
of the transmitted signal plus the Doppler spread of 
the reflected signal is the total observed Doppler 
spread. Thus, the Doppler spread due to the moon's 
libration is 

1 A/S=y c ^=[flin K2^-^a-2W-cos v(2l' L -Z La -Z Lb )}. 

Hence with the use of eqs (33) and (34) this may be 
expressed as 

Af d =2f c ^(smv-l Tfi --cosv.l TL ). (36) 

> c 

There is a small additional Doppler spread term 
which arises from the fact that the observer is 
approaching the general point L on the moon's sur- 
face, rather than its center. The original main 
Doppler shift, given by eq (21), was developed for 
j the center point of the moon. Actually it is the rate 
of change of D m in figure 10 which should be con- 
sidered, i.e., D m =D-\- small correction, and this 
correction will be different for the transmitting and 
receiving stations. Now D m ^D may be expressed 
&&'D m =~D')Jl — (k 2 Rm/Dm). The first term expansion 
of D m is sufficiently accurate for our purposes, so 
that with D substituted for D m in the expansion term 
the correction is given by 



D n 



1 • h 2 T? 2 

-v 2 u D2 



Thus the additional Doppler spread term is 

■Ml (37> 



a/I-i^F. 



where the subscripts a and b refer to the two stations. 
The total Doppler spread is now given by 

Af d =2f c — - [sin v> l T $— cos v- l TL ] 



+ l J ~l 



2 c'^ tm lDl + Dli (38) 
The maximum (or minimum) value of the Doppler 



spread occurs at the angle v given by 
p =arctan [— (l T p/l TL )]. 



(39) 



In order to obtain the maximum Doppler spread 
value, the positive value of (36) is chosen and added 
!o the absolute magnitude of (37). 

When the libration rates are not too small in mag- 
nitude (for example, not less than 10~ 7 rad/sec) a 
number of approximations can be made in the com- 
putations. In eqs (31) and (32), terms involving 
8 m and Di/DxX^l can be ignored. The term A/i' 
(eq (37)) also may be neglected. Check calculations 
show that in general a negligible effect is produced by 
these approximations. For small libration rates their 
effect needs to be included. 

However, it is to be noted that the calculations 
for the Doppler spread involve the daily values of 
selenographic latitude and longitude. The American 
Ephemeris and Nautical Almanac lists these to only 
three significant figures, so that at some point it is 
to be expected that lack of accuracy in these data 
will affect the computations. This appears to be 
the case at periods of very low libration rate. 

4.7. Relationship of "Fading Rate" and "Maximum 
Doppler Spread' 1 

The "fading rate" of the detected signal was de- 
termined as outlined in section 3, i.e., by counting 
the number of maximums in a given time interval. 
This obviously gives some measure of the form of 
the returned signal spectrum. It is apparent, phys- 
ically, that the total signal is due to reflections tak- 
ing place over much of the moon's surface, probably 
with a larger weighting factor attached to the cen- 
tral regions. If the spectrum is regarded as similar 
to that obtained when Gaussian noise is passed 
through a narrow Gaussian bandpass filter, then 
work by Rice [4] indicates that the expected number 
of maximums per second is 

N~2.52o- 

where a is the standard deviation of the Gaussian 
filter. Suppose the noise bandwidth, B c/s, is defined 
as that bandwidth containing half the total power 
passed by the filter, then 



or 



5~1.35 o- 



5^0.54 JVc/s. 



This gives a relationship between the bandwidth 
and fading rate that may be useful as applied to the 
lunar signal. The maximum Doppler spread (Af d ) max 
is a measure of the change in Doppler experi- 
enced at the limb of the moon due to libration. The 
maximum band limits of the signal, if its spreading 
were entirely due to libration effects, would be 2A/ d 
corresponding to the range of signal returns from 
opposite limbs of the moon. 
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The fading rate of the lunar signal was found to 
be proportional to the libration rate and hence to 
the value of (A/ d ) max . The best fit to the data 
was found to be given with k=Q.Q7 by the empirical 
formula 



N=k (Af d ) max = 0.Q7 (A/ d )„ 



(40) 



If our assumed relationship between fading rate 
and bandwidth is correct, then 



5=0.54 Xk (4/^=0.36 (A/,), 



(41) 



The bandwidth obtained is 0.36/2 = 0.18 of the 
maximum possible band limits due to the libration 
spreading. Since the change in Doppler Af d is pro- 
portional to the distance from the center of the 
moon's disk, then it is apparent that most of the 
signal energy is reflected from a portion of the moon 
about the center of the face with a fractional radius 
of only about 0.2. 

4.8. Actual Computations 

As results were required at intervals of a few 
minutes over a period of months, the computations 
were programmed for the IBM 704 and 709 Elec- 
tronic Data-Processing Machines. Much of the 
astronomical input data for the program was ob- 
tained on standard IBM cards from the U.S. Naval 
Observatory. Other input quantities were prepared 
from the American Ephemeris. The machines also 
were used for additional computations concerned 
with the antenna orientations. 

5. Results Obtained 

Figures 11 (a) to (g) show plots of some of the 
measured fading rates. The computed curves are 
based on a value of 0.67 in the empirical relation- 
ship of (40). The agreement in general is excellent, 
except at periods of very low fading rate of the order 
of 0.05 to 0.005 c/s. 

In general, the computed values for August 1957 
did not fall below approximately 0.1 c/s for either the 
approximate or the more accurate methods of compu- 
tation. The approximations involved in the method 
of computation are capable of allowing a better 
resolution than this, and it appears that the limiting 
factor probably is due to the 3-figure accuracy of the 
published daily values of selenographic latitude and 
longitude given in the American Ephemeris and 
Nautical Almanac. 

The value of k used, 0.67, when inserted in eq (41) 
for the "noise bandwidth" B, gives 

£=0.36 (A/,) max c/s. 

On this basis the major portion of the reflected signal 
is obtained from a central disk on the moon, the 
fractional radius of this disk being approximately 
0.2. A different approach is afforded by the results 
of experiments carried out at the Naval Research 
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Figure 11 a, b, c, d, e, f, g. Comparison of measured and 
computed fading rates. 
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Laboratory [5, 6], in which the interpretation of the 
form of reflected pulsed signals was that only a 
relatively small central portion of the moon's surface 
was effective as a reflector. 

Although not described here, it is of interest to 
note that the measured total signal levels showed no 
evidence of any slow fades. The Faraday rotation 
was observed, but the total received signal level, 
after correction for the varying earth-moon distance, 
was essentially constant. 

6. Conclusions 

Calculation of the rapid fading rate of moon-re- 
flected signals, in the manner described, has been 
shown to give results which agree very well with 
measured values. It appears as if most of the re- 
flected signal may be obtained from a central disk 
on the moon with a fractional radius of approximately 
0.2. 

The variation of the position of the region on the 
moon's surface which gives the maximum Doppler 
spread is of interest. Table 1 lists values of the 
angle v {) for a typical day, together with the angles 
of elevation and azimuth. The angle v depends 
upon the rates of change of the moon's selenographic 
latitude and longitude, so that the figures given in 
table 1 vary from day to day. However, the varia- 
tion is typical; values of approximately 200° to 300° 
during moonrise, an increase to values of 10° to 30° 
while going through transit, and then a further in- 
crease to 70° to 130° during moonset. It is apparent 
that whatever the value of the fractional radius of 
the moon is taken to be, a considerable portion of the 
moon's surface is covered during the course of a day. 
Since the agreement between the computed and 
measured fading rates is so striking, it would seem 
that, in its behavior as a reflector or scatterer at 
400 Mc/s, the moon's surface is reasonably uniform 
over quite a large region. 

Table 1. Variation of angle v with time for 21 August 1957 
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From the measured signal levels and a knowledge 
of the parameters of the system, the usual radar 
equation was applied to give an estimate of the 



cross section of the moon. This was 

o-=7X10 n m 2 . 

The projected area of the moon is approximately 
0.95X10 13 m 2 , so that the radar cross section, at 
412.85 Mc/s, represents 0.074 of the projected area. 
The manner in which the effective cross section of 
the moon may be proportioned between "gain" and 
"reflectivity" figures has been the subject of some dis- 
cussion[2, 7, 8]. If an element of surface is assumed 
to scatter in a fairly direction;) 1 manner, for example 
according to (cos d)'% where typical k values may 
be 30, 40, 50, etc., then most of the return is given 
by the ceii I nil portion of the moon, and the "gain" 
of the moon as a reflector, relative to its cross 
sectional area, approaches 4. With this sort of 
model, which is not critical at all, the power reflec- 
tivity must be of the order of 0.02. Values of 0.1 
have been quoted, which appear to be high, even 
allowing for a possible 3-db error in the measure- 
ments. 

The predictability of the rapid fading rate of a 
moon-reflected signal is of interest when considering 
the effects of the ionosphere upon transmission of a 
UHF signal. The signal is fairly well described 
when its average amplitude is known, the Faraday 
rotation allowed for, the Doppler shift computed, 
and the fading rate predicted to within a constant 
factor. Changes in the transmission medium, for 
example, auroral effects, then may be investigated 
by measuring their effects upon the known char- 
acteristics of the undisturbed signal. As mentioned 
in the introduction, such an experiment has been 
carried out between Alaska and Massachusetts. 
This experiment will be described in a forthcoming 
article. 

The authors thank the members of Lincoln 
Laboratory who contributed to this experiment. 
They are particularly indebted to S. C. Wang for 
his discussion of the libration calculations, and to 
J. C. James for his assistance in reading and correct- 
ing the manuscript. 
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